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' We argue that the van Dam-Veltman-Zakharov discontinuity arising in the hmit of the massive 

J> graviton through an expHcit Pauh-Fierz mass term could be absent in anti de Sitter space. This is 
possible if the graviton can acquire mass spontaneously from the higher curvature terms or/and the 
massless limit — > is attained faster than the cosmological constant A — > 0. We discuss the effects 
OO ■ of higher-curvature couplings and of an explicit cosmological term (A) on stability of such continuity 
and of massive excitations. 
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1 Introduction 



Has the graviton exactly zero-mass or perhaps a smah but non-zero mass? This is somehow related 
to whether or not the — > limit in a massive gravity corresponds to a strictly massless graviton, 
with M"^ being the mass of a spin-2 excitation. In the massive gravity theory defined by the Pauli- 
Fierz [|| term, one finds a discrete difference between the propagators for the massless graviton and 
that for the graviton with mass — > 0. Recently this has got a new impetus in an anti de Sitter 
{AdS) space j§, ^. The issue is subtle and appealing, and one has to resolve this by looking at the 
massive gravity theories defined by other than the Pauli-Fierz action. A key word of this approach 
is that van Dam-Veltman-Zakharov (vDVZ) discontinuity Q known to exist in Minkowski space may 
not survive in AdS space, in particular, if a smooth limit M^/A — > exists. Kogan et.al. in and 
Porrati in used the spin-2 Pauli-Fierz mass term to demonstrate the absence of vDVZ discontinuity 
in AdS space at tree level, and further realization of this non-perturbative continuity for the spin 3/2 
field oi J\f = lAdS supergravity has been given in ^ Absence of vDVZ discontinuity in de Sitter 
(dS) space was already shown by Higuchi ||^], but as the result was restricted there to de Sitter space 
and hence may have less importance, for there is no unitary spin-2 representation in the mass range 
< < 2A/3. A smooth limit M^/A has been observed in Q for the four-dimensional theory 
of massive AdS gravity by embedding it into five dimensional AdS brane-world models. 

In order to achieve a smooth limit M^/A 0, perhaps one should find a mass term of the order 
A^, which could arise from the terms that are quadratic in the curvature. Naively, one finds reasonable 
to introduce the higher curvature terms in a ghost free combination into the starting action. A concrete 
example is the Gauss-Bonnet (GB) term j|9|, |lO|, which is a ghost-free combination (see Ref. [ |Tl| ] for 
the absence of ghost state in AdS^ based a single brane model). But it is obvious that the GB term 
does not generate a mass term for graviton in four-dimensional flat or curved spaces. Eventually, one 
could attempt the problem by introducing the higher-curvature terms in generic form. In this paper we 
explore this possibility and find a smooth limit for M^/A — > 0, where the could be quadratic one of 
A. It is possible that a very light graviton in the presence of non-zero cosmological constant does not 
show a significant deviation from the observational effects that of massless gravitons at large distance 
scales. For example, the prediction of Einstein's gravity for the bending of light by the Sun agrees 
to the observation with less than 1% difference, a limit M^/A < 10~^, therefore, must be satisfied 
by the graviton mass M^, and hence a small graviton mass may not be in contradiction with the 
observation |12, 13 1. 

An undeniable fact is that the starting gravitational theory is Einstein gravity, thus any natural 
modification on it should arise from the higher order corrections in the curvature. Among the different 
motivations in higher derivative theory, the special one here is attributed to an expectation that mass 
term arises from terms that are quadratic in the curvature, so that a smooth M^/A limit can 
exist. It is really astonishing that four-derivatives gravitation in a curved background that we study 
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here is a crucial tool to show the absence of van Dam-Veltman-Zakharov discontinuity in an AdS space. 
Without any ambiguity we understand that the higher derivative theory with terms upto quadratic in 
the curvatures effectively contains, which also we explicitly demonstrate in this paper, in addition to 
the usual massless graviton excitation, a massive spin-two ghost and a massive scalar. Nonetheless, by 
judiciously choosing the higher derivative coupling /? (= /? + 47) zero, one can always recover a theory 
which is ghost free. The most radical result is that one recovers the correct expressions for the massless 
and the massive graviton propagators by taking simultaneously the limits /? — > 0, A ^ but A >> /?, 
and A — > 0, /3 — > but A « (3. We also present several novel results in a flat spacetime background. 

The outline of the paper is as follows. In Section 2 we give some basic features of the fourth-order 
gravity in a curved background showing that one generates a massless spin-2 graviton and extra massive 
degrees of freedom with the higher derivatives. In Section 3 we evaluate the one-particle exchange 
amplitude between two conserved sources and demonstrate the absence of vDVZ discontinuity in the 
A — > limit. In Section 4 we briefly discuss the issue of the general covariance with a (non-)trivial 
A. We also briefly describe the spin and field content of the linearized four-derivative equations, and 
the counting of degrees of freedom by introducing gauge fixing term. Section 5 consists of conclusions. 
Finally, Appendix A presents the basic variational formulae for different curvature terms. 

2 Effective action and Linear equations 

We start with the Einstein-Hilbert action supplemented with the quadratic curvature terms: 

Sg = I d^x^{ ^R-2A + aR^ + /3R% + jR^p^} + Sgf . (1) 

Here Sgf is the gauge fixing action. Here K^a, have dimensions of inverse mass squared, and 

K?K has dimension M^. For simplicity, we set = I^ttG^ = 1, unless otherwise stated. We may keep 
A arbitrary, but most of the results bear meaningful interpretations for A < or A = 0. It is known 
that the theory defined by (||) explains, in addition to the usual massless graviton, two new fields: a 
massive spin-2 field with five degrees of freedom (DOF), and a massive scalar with one DOF |T^ , |T5| . 
The equations of motion derived from (|l]) read as in units = 1, 

{R^.u - \g^..R^ + Ag^, + 2aR (^R^, - + 2/3 {r^p.^R"" - -^g^uRp^RP""^ + 

+ 27 (^Rf^paSRu — ^ QiivRpa&xR'"^^^ — 2-R^ ^Rv\ + 2Rp,XypR'^P 

-{2a + [5 + 27) (V^V, - gpuV^) R + {(3 + 47) (^Rp, - 1 gp,R^ = (2) 

In the linear approximations, the curvature terms like RRpi,, Rppu^R^" affect the gravitational exci- 
tation spectrum near flat space. That is, in a flat spacetime background (A = 0), second, third and 
fourth bracket terms in (|2p do not contribute to the equations of motion linear in fluctuation, but they 
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do contribute in the curved backgrounds. The AdS (or dS) space solution is expressed in terms of the 
background metric (g^u) as 

Ri^iupa = -^{g^lp9va — 9na9up), Rfiu = Kg^y, R = 4A. (3) 

Here A = —3//^ is defined, where / is the AdS4, length scale. It is trivial to check that the background 
solution (|3|) also solves the zeroth order equations of motion (e.o.m.) of higher derivative gravity 
defined by (|l|). One must replace gfj_i, in (Q) by g^^- Of course, the covariant derivatives taken about 
the constant (negative) curvature backgrounds are vanishing. A field satisfying R^^, = Agf^i^, R = AA 
always solves (§). The sum of the first two terms (i.e. G^i, + A^^j^) in (§) gives zero. Also note 
that second, third, and fourth round bracket term in (§) taken about the background (^), is each 
separately vanishing, hence one finds the standard e.o.m. of massless gravity for any value of A. This 
is completely a consistent argument, because unless we perturb the background, we do not except any 
massive (tensor or scalar) modes to be present. 

In order to study the propagation of the metric and the complete mass spectrum in four-dimensions, 
we introduce the perturbation around the background space gfj^iy = g^u+h^jy. Then to the order linear in 
h^y, neither the derivatives terms, nor the bracket terms in (^) are vanishing in the curved backgrounds. 
Note that the source term for the metric fluctuation h^i, can then be defined by 



Sjvi = —2 



j d^x^gh^.T^'^ . (4) 



We can further make a simplifying assumption that Tfj_y is co-variantly conserved with respect to 
the background metric, Vf^T^'^ = 0. For the study of massive gravity, a common practice in the 
literature [1-7,17] is that one adds the following Pauli-Fierz action to the linearized action for h^i, 
derived from (|I|) (but setting a = /3 = 7 = 0) 



SpF = Mlp / d^x 



Kuh^'' - [hx^f] , (5) 

where Mp^, is a pure spin-2 mass term. A clear motivation for introducing this term is that it is a 
ghost-free term for a free, massive spin-2 propagation at the linearized level. However, as the mass is 
acquired in an explicit way, one cannot, in general, expect a smooth limit at the quantum level |17|. 



Remarkably, we find that (^) is not essential to prove the absence of discontinuity in graviton mass 
in an AdS background. Beside, we have no good reason yet why a spin-2 excitation should have an 
arbitrary real mass in the pure four-dimensional gravity. Indeed, in a pure 4d theory, addition of the 
leading order quadratic curvature terms to the Einstein-Hilbert action (the total action then explains 
a massive or multi mass model of gravity) is more logical rather than an addition of the explicit Pauli- 
Fierz (PF) mass term. Besides, a purely four-dimensional theory of massive gravity with the PF term 
is not well defined in four-dimensions [|l^], because an arbitrary spin-2 field does not guarantee that this 
field will couple to a conserved source. The origin of this term is expected to arise from some massive 
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graviton sector of the usual Kaluza-Klein reduction (see for example Ref. [18| for a close realization). 
We simply the drop the Pauli-Fierz action in the present analysis. 

In order to obtain linear equations of motion for h^i, we use the following relations 

25Rf,, = V- VmV./i + 2V(^VX)p, (7) 
5R = -R^^h^'f -V^h + V^Vph'^f , (8) 

(2) 

where the Lichnerowicz operator A]^ acting over the field h^^u is given by 

Ai^^ V = -V^ V - ^Rpf.auhP" + 2Rp^phP ^) = -V2 V + ^ ( V - \ 9nvh^ , (9) 
and the operators A^'s further obey the following properties [19| 

Af = -g^uV^h , A^^Vm = (-V2 + A) . (10) 

Here A^^^ and A^^ are the Lichnerowicz (differential) operators acting respectively on the spin-1 and 
spin-0 fields. Note that A^-* , A^^^ and A^^ commute with contraction and covariant differentiation 
about the background. We can easily show that the local curvature squared terms in the GB combi- 
nation does not generate mass for the linear perturbations. This is because in ordinary two derivative 
gravity supplemented by the GB term, each of the four gauge-group local parameters of Diff-invariance 
accounts for killing of two degrees of freedom (DOF), leaving just two DOF of the massless graviton 
out of ten DOF of h^^u. A general feature of gravitational action with arbitrary a, /3,7 is, however, 
that the graviton excitations near a flat space could behave as ghosts. This means that the set of all 
states would not form a Hilbert space with a positive definite metric. However, it is worth noting that 
the theory defined by (|l|) is ghost free when a relation /3 = —47 holds. 

To the first order in h^i, and up to the terms that vanish on the Einstein backgrounds (|^ (see 
Appendix A for the details), the linear equations of motion in D = 4 take the following form 

(1 + ai)A^^^ V + 2(1 + a2)V(^V^/i^)^ - (1 + a3)V(^V^)/i + 2A/i^^ 
-4(1 + a3)A( V - -^9f.uh) + (1 + a4)g^,uy^h - (1 + ar,)g^^V ph^P 

+ 2(2a + (3){VpV, - gp,V^) (V^/i - VaV^/i^'') = AT^, , (11) 

where is taken to be covariant with V^gxp = 0, Q; = a — 7, /? = /? + 47, and the coefficients Oj are 
defined by 

ai = 8Aa + /3fv2- 02 = 8Aa + /^fv^ + ^ 
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,2 . 2A^ 



as =4Aa + /3(V^ + — j . (12) 



The trace of Eq. (O) takes the following form 



2-4/3 ^V^ + t) ~ + 



- 4/3V^ - 12aV2 VVhf,^ + 2A/i = 4r . (13) 
An obvious consequence of fourth-order derivatives of the metric in field equations is that there are 



six more degrees of freedom in the theory other than the two massless graviton states [14]. One may 
also define a condition to constrain the number of propagating degrees of freedom. A physical choice 
is V^V'^/i^jy = V^/i, which indeed satisfies the field decompositions given in the Ref. [0]. It is worth 
noting that in the massive gravity theory given by the Pauli-Fierz(PF) action, this constraint naturally 
appears from the field equations j^, which simply follows by taking a double divergence of the linearized 
equation containing Mpp{h^^ — g^uh). Since the covariant derivative of the Einstein-Hilbert term is 
zero, one finds Mlp{yh^„ - Vyh) = 0. This further implies that V^V/i^,, = V^/i for on-shell h^y, 
and as we see below this is more useful to compute the one-particle amplitude. 
Then from the Eq. (|l^), we obtain 



(37W2 _ 2A)/i = -4T, (14) 
where we have defined M'^ = 2K((3 + 6d) VV3. Meanwhile, if we choose /3 + 3a = in (|l|), we find 



- 2V2/i + 2V''V^ V - 2A/i + ^M^h = -AT . (15) 

By replacing M."^ with Mpp, one observes that the Eqs. (|l^ ) and ([l5| ) become identical with the 
Eqs.(13) and (11) of the Ref. Nevertheless, there is a clear difference in their origins, and Ai^ 
defined above is not what a PF mass term is. Mpp is a pure mass term put by hand, while A^^ is 
just like a mass parameter, whose value is determined only when one defines the pole V^. When we 
express the curvature terms in the GB combination, we have d = /3 = 0, and hence = 0. This 
is expected because the GB combination is a topological in four-dimensions and does not change the 
Einstein's gravity. 



To get something in concrete form from the Eqs. ( [ll| ) and (13), we need to deploy the following 
relations 

(v,Vx + I g,x^ / = V, (^V^Va + | g^x^ f (16) 
^2 ^ 2A \ „ ^ ^ 8A 



(^V^V, + — g^,j f = V^V, ^V' + —jf (17) 
V'^ (V^/. + V,/^) = -Al^fu + 2AU + ■ (18) 

Apply to Eq. (|Tl|). Using the above relations we find 

A/3 (V^ - llA) V - A/3V^V'^V" V = • (19) 
For instance, if we impose V^V'^/i^o- = 0, then Eq. ([T9| ) implies that 

A/3 (V^ - llA) V'' V = • (20) 
5 



Here if A/3 (V^ — llA) is defined by some mass parameter m^, the Eq. (pO|) reveals that only 10 — 4 = 6 
(Not 5!) components of h^y propagate independently. Five of these correspond to a massive spin-2 
particle and the sixth corresponds to a massive scalar. However, with V^V'^h^^ = V^/i, we can further 
reduce the number to the five DOF of a massive spin-2 field, which corresponds to a ghost-like particle. 
But we can make theory free from ghost by choosing /? = 0, which is the choice we would like to have 
in most of the discussion. 

Finally, compute the double divergence of Eq. (|ll|) (or simply apply V'^ to Eq. (p!9|)) and arrive to 



- lOA^/^V'^V'' V = • (21) 

As we see below, for the absence of an unphysical pole, we must take |A/3| = 0, so the possible choices 
are A = 0, /3 7^ 0; A 7^ 0, /? = 0; A = /3 = 0. The choice /3 = is more preferred because it makes 
the theory free from the massive spin-2 ghost, and also makes the residue at an unphysical pole zero. 
Otherwise the equality V'^V^h^u = should hold. 

3 One-particle amplitude 

Now we decompose h^u as follows 

Ku = hl^ + Vi^^V^) + V^Sv(^ + 9^ui^- (22) 
This, with the conditions ^V^ = and V'^/i^J = = K^'^ , further implies that 

v2/i = vV + 4VV, V^V'^V = V^(/> + AV2(/> + VV- (23) 

The one-particle amplitude is most easily computed by choosing V^V^hxp = V^/i, thus from (^) we 
find SV^V = AX/'^cf). From (^), we then easily read off /i = (3V^ + AA)ijj/ A. Use this value of h in 
Eq. (0) and arrive to 

(3V^ + 4A) = -4A (2A (p + 6d) - 2a) ~ T (24) 



The transverse traceless (TT) part of Eq. (11) is 

\l + ai)Af - 2A(1 + 2as)) hJ^J = 4rJJ , (25) 
where the TT-component of T^,y can be expressed in the following form 



= V - ^ 9f.uT + ^ (V^V, + gp,A/3) (V^ + 4A/3) ' T . (26) 



Then the one-particle exchange amplitude between two co-variantly conserved sources T'^^, and T^,y can 
be simply written as 
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2T 



/flU 



1 



^l + ai)A^^^ -2A(l + 2a3; 

1 



. T - - T' ^ 

^"^ 3 -(l + ai)V2-2A(l + 2a3) 



T 



(-(l + ai)V2-2A(l + 2a3)) (V2 + 4A/3) 



2A_, 1 
T 

3 2A + 6d)V2 - (V2 + 4A/3) 



T. 



(27) 



Obviously, there are four propagators that may have poles in the amplitude Eq. (^), and we have to 
find the residue at each pole. We now proceed to study each and every poles. 

(A) Pole at V2 = -4A/3 

The residue at the unphysical pole = — 4A/3 is given by 



2A\ -3 

"9"; 2A + 8A2(2a-/3) 



2A 
~3~ 



-6A-8A2(6a + /3) 
16A/3 



1 + 4A(2a - (3) 3 + 4A(6d + (3) 



(28) 



Hence the residue can be made zero only by taking A = or /3 = 0, these are the choices we would 
like to have henceforth. 



(B) Pole at V' 



2A(l+2a3) 
(1+ai) 



As stated above, for a meaningful result, the residue at the unphysical pole should to zero, which is 
achieved by taking A = or /? = 0. The propagator structure, however, depends upon the fact that 
out of (3 and A which one vanishes first. Thus we find two possibilities. 

(I) |A| »/3 

By setting (3 = first, one kills the massive spin-2 degrees of freedom, and hence the usual massless 
spin-2 graviton propagator can be recovered. By substituting the values of ai and as from (p!2[), we 
obtain 



2A(1 + 03) 
(1 + ai) 



-2A 



l + 8Ad + 2/3 (V2-2A/3) 



1 + 8Ad + /3 (V2 - 8A/3) 



-2Ae// + M2_^^ = -2A, 



where 



A, 



eff 



A_ 

1 + 8Ad 



Mi 



2M^ 
1 + 8Ad 



8dA^ > 0, (d < 0, |Ad| << 1) . 



(29) 



(30) 
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The residue at = is 

r 9 9A /— ?<M 1 1 

(31) 



2 _^ 2A /-3\ 



3 9 V2A/J l + 8Aa 1 + 8Aa 
In the hmit A — , i.e. — >■ 0, the residue is always —1, so that 

^ = iTo (OT) °" - ^' °" ^1 • 

where □ = -V^ - 2Ae// + M^^^. It might be relevant to compare this to the physical pole at 
= M'^pp - 2A when one introduces an explicit Pauli-Fierz mass term instead of higher derivative 
terms. When d ^ 0, Ag// A and Afgjj 0. Further, if we set A = 0, one finds the well known one- 
particle exchange amplitude for a massless graviton in the pure Einstein's theory. This may be enough 
to understand the absence of vDVZ discontinuity in an AdS space. As expected, for /3 = 0, we find a 
massless graviton of the Einstein-Hilbert theory, but the 4d Newton constant can be renornializcd due 
to the non-trivial a. Since all the extra degrees of freedom decouple for /3 = 0, there is no any ghost 
state in the theory, and also that the residue at an unphysical pole mentioned above is vanishing. 

(II) |A| « ~^ 

For general a. and /3, the two physical poles at = imply the following quadratic equation 

^V^ + (^1 + 8Aa + j + 2A + 8A« - ^ j = . (33) 
Clearly, there are two poles corresponding to positive and negative roots of (with |dA| << 1 and 

\m « 1) 

8A^ / \ 

Vl = -2A+ — (f5-6a) =-2A + Ml, (34) 



4A (1 + 8Aa) _ 4A 



n = -^-'-^^^^-^-Ml, (35) 



If A goes to zero first, we have Mi = and M2 = 1/(3- If we further send /3 — 0, the massive spin-2 
field decouples from the theory. Thus, at the negative root pole V?. = — 2A -|- Mi, the residue is 

"I 7 7^ \Y ' 

^ 3fl-4Af/?-6aJJ 

which is always —1 in the limit A — > 0, so that 

A = Gn [2r""' T^u - T' T] , (37) 

where V = j3 {v"^ + l/jS^ V^. While at the positive root pole V+ = - Ml, the residue is 

_2 2A p 

S + T' (l + 8Ad)' ^^^^ 



which is always —2/3 in the hmit A ^ 0, so that 



A = Gn 



- T'V^^ T 
3 



(39) 



The results and p9| ) are the well known the vDVZ amplitudes for massless and massive spin-2 
fields in the flat space-time backgrounds. In particular, for /? 7^ 0, the lowest mass of a localized 
graviton in an AdS^ space is quadratic one of A, i.e., SM"^ oc A^. The amplitude ( [39| ) remains the same 
even when we take the limit /3 = 0, except that V = V^. So a clear message with the higher derivatives 
is that, even though the propagator is found to have a single pole in the limit A — > 0, /? — > 0, (i.e., 
V = 1/V^), there are two different propagators - one is massless (Eq. (p7|)) and another is massive 
(Eq. (p9D). This has been possible due to two different root solutions for V^. Obviously there is no 
spin-2 ghost with /3 = and also the residue at = — 4A/3 becomes zero. 



(C). Pole at {(5 + 6a) -1 = 



In fact the amplitude (27) is singular at = (/3-|-6a) ^. The amplitude is sing ular at M^p = 2 A/3 m 
the massive gravity defined by the Pauli-Fierz action This allows us to define A^^ = — 2A(/3-|-6a)V2, 
and the above singularity is similar to the singularity at M'pp = 2A/3. However, the amplitude is 
smooth in the limit < A^^ < 2A/3, and a smooth limit TW^/A ^ is possible to be achieved since 
~ A^. A more remarkable point is that when we set d = /3 = 0, so that the curvature squared 
terms are in Gauss-Bonnet combination, neither there is any unphysical pole nor there is any ghost. 
Nonetheless we recover correct amplitudes for both the massless and massive gravitons by taking, 
respectively, /3 = 0, A — > 0, and A ^ /3 — > limits. This is a quite remarkable result. 

3.1 Flat space limit A = 

We find useful to discuss some important features of the four-derivatives gravitation by taking the flat 
space (A = 0) limit. The non-relativistic potential for arbitrary a, /3,7 takes the following form (taking 



(40) 



167rG4) p4|, pof (see references in pG] for earlier work) 



U{r) = GiM^ 



where 



ml 



1 

r 

1 



2 g-mor ^ g-m.2r 

3 r 3 r 



1 



TTT-n 



(41) 



2^2 1^/3 + 3aj ' k'^P 

and, mo, rn2 are the mass terms for spin-0 and spin-2 particles. These massive modes would rise to 
give Yukawa-type interactions. With the choice P + 3a = , \(3A\ « 1, we find 



4A2/3 16A3/32 



+ 



2A 16A2/3 



(42) 
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In the limit /3 — > 0, = 2 A/3 and Ai^ = 0. Clearly, when Ai'^ represents mass of the spin-2 
graviton of pure Einstein's theory, may characterize a mass term for the graviton, but this also 
vanishes for A = 0. In the limit A — > 0, though both and approach zero, A4?_ reaches to 
zero faster than A^^. For /3 7^ 0, the massive spin-2 field can behave as ghost-like. Another way of 
seeing the presence of ghost-behavior is that one of the correction terms to the Newtonian potential 
in ( |40| ) has a coefficient 4/3 with wrong (positive) sign, and an ultimate way to make the theory free 
of ghost at the tree level is to set /? = 0. It is remarkable to note, which is also recently reviewed 



in 21 1, that for finite mo, m2 and at large distances, the potential (|4y) reduces to Newtonian limit, 
and it is finite for r ^ 0. It was known that there is no vDVZ discontinuity even with the Pauli-Fierz 
term when one sums up all the tree-level graphs i.e., in the full non-linear analysis. Thus, with a 
non-zero cosmological constant, regardless one introduces Pauli-Fierz term or higher derivatives, there 
is a continuous limit for ^ in the lowest order tree level approximation, and we argue that this 
is a property of any space with non-zero curvature. 
From Eq. (|l]) one finds 

ml + ml= ^ + , >0. (43) 
2/3 (f3 + 3a 



When (3 + 6a = 0, one has niQ = —m\ = 1//3. On the other hand, in the same limit, one has 
= = = A'P. Since = limit corresponds to the spin-2 massless graviton, we can 
therefore take the limits a — > , /5 ^ (i.e. mo oc 00, m2 oc 00) to arrive at the standard Newtonian 
limit for non-relativistic potential. 

A comment is in order. As stated above, for a non-zero (3, there may arise the problem of negative 
probabilities for processes involving an odd number of massive spin-2 quanta and hence violation 
of causality. However, since the parameters a and (3 are small enough to make the massive fields, 
the breakdown of causality might occur only on a microscopic scale or near a Planck scale. So a clear 
message is that if the ghost are gentle and harmless creature, with whom we can live quite comfortably, 
a choice /? 7^ 0, A = is possible, otherwise we must have /? = 0. As we have seen above, ghost are 
gregarious which demand the company of real particles, and they appear only at high energy, so their 
effects will be insignificant in normal particle scattering, which will appear to be unitary to a high 



degree of accuracy [23|. 



4 Gauge covariance and gauge fixing 

Any viable gauge in a general four-derivative gravity theory can involve third derivatives of the po- 
tentials (or metric fluctuations), unlike the de Bonder gauge in the Einstein-Hilbert theory, which 
involves only first derivatives of the potential. In a flat space background A = 0, the linearized field 
equation (|l^) is manifestly invariant under the coordinate transformation — > x'^ = x^ -\-^^, where 
^^{x) is infinitesimal vector field. For A 7^ 0, however, the metric fluctuations h^^, are not, by them- 
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selves, invariant under infinitesimal transformation. Here we shall make a cursory inspection of gauge 
invariance for the A 7^ and A = cases. 

4.1 A non- vanishing A 

To justify the above mentioned argument, let's write the action in the following form (i.e., setting 
= 1, /? = 0, 7 = 0, so that /? = and d = a) 

h = j d^x^/^ {R-2A + aR^) - 2 j d'^x^/^h^^T^'' . (44) 

For the background solutions (^, the equations linear in h^y read 

^i?(V) + 2AV - AK^a (4V - g^^/i) + SAaA^^h^^ 

+^Aag^y - VxVph^") + 8Aa (2V(^V^/i^)a - V^V^/i) 

-4a (V^V, - g^,V2) (-V^h + V^Vph^f - Ah) = AT^, , (45) 



where is the linearized Einstein term given by 

V + 2V(^V^/i,)^ - V(^V,)/i - A(4V - g^,h) + g^, (v^h - VxVph^") (46) 



As is known, a = explains only the spin-2 massless graviton, thus the Eq. (^5|) is clearly invariant 
under the gauge transformation h^i, h^y + V [pV^y One can use this invariance to set 

V'^V = ^V,/i, (47) 
which removes four degrees of freedom. Set T^^ = and a = in Eq. (p5|), find its trace, multiply it by 



gpu/2, and subtract the resultant expression from the Eq. (45) itself. Then one can use the gauge (| 
in order to arrive at 

The residual gauge invariance is then used to fix the above gauge, which is solved by 

V^V(^K) - \ V.V'^V; = i (V2 + A) K = . (49) 

The Eq. (|8|) then implies 

(a« - 2a) Vy) = -V(^ (V2 + A) F,) = . (50) 



Since this equation is exactly solved by the gauge fixing condition ([49D, there remains only 2 (10 — 4 — 4) 
physical propagating degrees of freedom 1^], which correspond to a massless spin-2 graviton. 
Now for arbitrary q, using the same gauge transformations as above, we find 

(1 + 8Aa)V(^ (aJ;^ - 2A) V,) + 2aV^V, (V^ - A) V 

+ ag^y (V^ + 2A) (V^ + 2A) Va^^ = . (51) 
11 



The gauge- fixing condition does not automatically solve this equation unless a = or VxV'^ = 0. 
One finds then three^ (= 10 — 4 — 3 , not just two!) physical degrees of freedom in the theory defined 



by (44) - two of them correspond to a massless graviton and the rest one to a massive scalar excitation. 



B. Whitt in |25] showed that the fourth-order gravity theory defined by (44) is conformally equivalent 
to Einstein's gravity with a massive scalar field, thus there is an extra scalar DOF in the theory defined 
by (0). 

Further we consider the action by adding to Einstein-Hilbert action a term quadratic in Ricci tensor 

h = J (fx^ {R-2K + (3R^uR^n - 2 J d^x^ h^.T^''' , (52) 

The linearized equations for this action take the following form, up to the terms that vanish on Einstein 
backgrounds (^), 

(l - (3A^^^ + 2A/3) (V(^ (a^^) - 2a) )) - /3V^V, (aJ^^ - 2a) VxV^ 
-4A2/3V(,K) + f ) ( aSP) - 2a) V,V' - ^ g.^V.V' = . (53) 

Obviously, the gauge fixing condition (|4^) does not solve this equation. What this implicitly implies 
is that the theory defined by (|^) has massive modes, other than the usual massless spin-2 graviton, 
i.e., not all the components of /i^jy can be fixed by a single gauge. It is known that the theory defined 



by (^) has a massive spin-2 mode other than the usual massless graviton flj, |28 |. 
4.2 The vanishing A 

In a flat space background, one may separate the fluctuations h^^ into massless spin-2 field0(/>^j,y(= /i^^"*) 
and massive fields T,^^. The field can be broken up into a pure spin-2 field u^i/ and a massive 
scalar ijj. Being more precise 

h^,u = 4>^lv + '^^lv , '^^lv = cf^^ + rj^yip + 2m~'^d^duTp ■ (54) 

However, since we also have the massless field (pf^u, the gradient term above may be dropped, for it can 
be absorbed into a gauge transformation of Not surprisingly, the gauge invariance of the original 
action is characterized by a gauge invariance for the (pf^i, field alone. 

With A = 0, the linear equations derived in (|l^) take the following form 

1 + (^SR^, - ^ ri^,6R^ - (/3 + 2a) {d^d, - ri^.a) 6R = 2T^, , (55) 



where 



1 



n = V^iud'^d" , h^u = hpy- -ri^^h. (56) 

^Notice that there should be an extra constraint, so that only 3 out of the 4 gauge variables are free. 
■^Note that h']^^ represents the massless spin-2 field (graviton) of the Einstein's theory when a = /3 = 0. 
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The trace of Eq. (55) gives 



1 



l + pa]6R- ((3 + 2a] a6R 



T. 



Then from the Eqs. ( |55[ ) and (p7D, by using (p6[), we arrive to 



1 



where 

= (i + - (/3 + 2d) 5R . 

As expected, for d = /3 = 0, the field equations (|58|) and ([57|) in the gauge d^^h^p 

1 



yield 



□ 



-4 T, 



where hli^' is represented by (p^,^. Notice that under the coordinate transformation 
the variables r^(x) transform as 



Tf,ix) — > r' (x) = r^(a 



(57) 



(58) 



(59) 



(60) 



(61) 



Now the gauge must be fixed by choosing some arbitrary relations which determine the in terms of 
arbitrary variables and their derivatives. Practically, one can also choose a gauge that determines 
in terms of the given Cauchy data on the initial hypersurface. A simple viable gauge is r^(x) = , 



which is known as "Teyssandier" gauge [26|, and this gauge choice was justified in Ref. [24| by choosing 
7 = 0. Obviously, for d = /? = 0, Teyssandier gauge reduces to de Bonder gauge d^h^p = 0. With 



= 0, the Eq. (58) can be written as 



4 T, 



(62) 



where m| = 1//? has been used, and the massive spin-2 field cr^jy is defined by 

_ -nhpu - {l/2,)rjpu5R 



(63) 



Indeed, consistency of Eqs. (|62D and (|63| ) with T^ = 0, latter as a gauge condition, also requires that 
d^d^apy — □cj = 0, which is the reminiscent of the constraint variables V'^V'^hpi, — V^/i = in the full 
theory defined by (Pl). Eq. (^) may be written as, using Eq. (|63|) , 



□a, 



^hpiy - r]i^,y6R = 4 (Xpi, - ^ r]p„T 



However, in the gauge Tp = 0, the trace part of (|58|) yields 



1 



l + /3a]ah 



l + (3a]6R 



■T, 



(64) 



(65) 
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while the gauge r„ = itself implies 



1 + pa)ah+ (1 + pa) 5R - {0 + 2a)a) 5R = 0. 



1 

2 

Hence the Eqs.(|65D and ( [6^ ) would rise to give 

l-2('/3 + 3d) □ 



6R = -2T . 



(66) 



(67) 



Substitute this value of 6R into the Eq. (64) and arrive to 

1 



3mQ 



r|^Jip 6R 



-4 [T^v-i^ ri^^T ) . 



where as defined previously rri^ = — 1/(2(/3 + 3d;)). Evidently, Eqs. (|60|) and (|68|) reveal that 



where the massive scalar field V' = 5R/{'im^) satisfies the following equation (from Eq. (|67D) 

(□+m2)V' = -^r. 



(69) 



(70) 



Thus, in the momentum space, the interaction amplitude between two conserved sources T'^^ and T^y 
is given by 



1 



^rp rpltlf _ 1 rpXrpIP^ ^rp rpffJ.l' _ 1 rpXrp/p^ rpXrpfP 



+ 



6 (p^ + ttIq) 



(71) 



Thus to make the theory free from a massive spin 2 ghost state, we set /3 = 0, so that the second term 
above is formally absent. 

4.3 Gauge-fixing term 

The gauge invariance of ( p5[ ) is made manifest by writing the action in the following bilinear form 

1 



S 



ct'x 



where the differential operator kernel jQ^"" for the diff-invariant part is 

= 0(1 + /?□) - 2a (l - 20 + 3d)n) Pl%, , 
where the transverse spin-2 and spin-0 (scalar) projectors are defined by 

^ fiuG pa I Vpi' G pu — 



p(2) -Ifl lifl e _\q e p{s) _ 1 



pi/pa ^ 2 "PP '^^ ^ 2 '^P 2 "pVJpa 1 



□ 



(72) 



(73) 



(74) 



Here P^^) ^nd P^'^) are symmetric under /i z/, /) <-> cr, /xz/ <-> po". As we illustrated in the last 
subsection, in the absence of gauge fixing term, K'^"-'" explains eight physical propagating degrees of 
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freedom in the spin subspace 2 © 5, which correspond to the massless graviton {(j)fj,u), the massive 
spin-2 field (cr^i/) and a massive scalar (ip). Indeed, we can make the four-derivative gravitational 
action complete by including to (72) the following gauge fixing term [27| 

Sgf= [ d'^xx^WG^uX'^ih), (75) 



where 

X^(/i) =9,/i^'^-AiV^/i, = A2V''Vp5M. + YV(^V,) + A4ff^., (76) 

with Ai,---,A4 being the gauge fixing parameters. The gauge fixing contribution in terms of the 
differential operator kernel is 

IC9f = -6X1 P^^^ - □ (A4 + OAs) P(^) - 2(1 - Ai) f (1 - Ai)nA* p(^) - AiQA* P^^^A (77) 



where A* = A4 + (A2 + As)^, and the spin-1 operator P^^^ and the transfer operators p(^), p{SW) ^^^.^ 
The total contribution to the four-derivative differential operator is given by the sum /C = IQ^"" -\-K,3f _ 



= - e^pu,„ + - e^^uj^p , P(^) = io^^iUp^ , = Op^iOp^ + Op^LOp, . (78) 



One can redefine the field hp^, by the following transformation [27| 

/C — >jC = AICA, (79) 

so that the operator ^(Ai) defined by 

^(Ai) = P(2) + p(i) + 2 ^(^) ^ 2 (1-Ai) ^(5,^) ^g^^ 

3 9 Ai 

is invertible, and becomes a numerical matrix for Ai = —2. While, the new field hpi, is transformed as 

hpu = iA-XZhp<. (81) 
Finally, the quartic propagator obtained by inverting the projectors take the following form p^ ] 

p(M = I p(2) _ ^ M 

(l + /3n)n 8 (i_Ai)2('i-2(/3 + 3a)n)n 



^ -piS). (82) 



(A4 + A2n)n 8 _ x^)^ (^x^ + (x^ + x^)n 



This can be split into two new fields h and vr as 

^'^^^ " 8 (i-Ai)2n^ A4n 8 (i-Ai)4A4n ' ^ ^ 



^ (D + mi) + 8 (1-Ai)2 2; 



+ . , ^ P(^) + ^ — ^ ^ P(^) . (84) 



A4(n + m2) 8 {l-Xir x,{a + m^) 



*As several authors have studied the fourth-order gravity action with 7 = supplemented by gauge fixing term (see 
Ref. for review) , we will be brief here. A theory defined with 7 7^ is more physical in a sense that one obtains ghost 
free action with the choice /3 + 47 = 0, but in the formal case one must set /3 = 0. 
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where 



mo = — , m,,, = = , mi = — , m = -— . (85) 

^ (3 2(/? + 3d) ^ Aa' ' (A2 + A3) ^ ^ 



The sign in front of each spin operator in Eqs. ( p3| ) and (84) is important. Obviously, in the massless h 
sector Eq. (p^), there are just two physical propagating degrees of freedom in P^^^, which correspond to 
the massless graviton, while other three ghost modes (if A4 > 0) are just gauge degrees. On the other 
hand, the massive vr sector (Eq. (|8^)) contains a massive spin-2 field of mass with five degrees of 
freedom and one physical scalar with mass equal to mg defined previously, and rest are just the gauge 
degrees of freedom, but they are not necessarily ghost excitations, it depends upon the sign of gauge 
parameter A4. If A4 < (provided also that A2 < and IA2I > IA3I), all gauge dependent massive fields 
in vr sector become real gauge dependent fields, otherwise they show ghost behavior. When we set 
/3 = and choose a gauge such that Ai = 0, we will be left only with a massless graviton as a physical 
spin 2 excitation of the Einstein's gravity. 

It is essential to quantize the theory in order to extract any concrete ideas about the vDVZ dis- 
continuity at the quantum level. Obviously, one must then add the compensating higher derivative 
Faddeev-Popov Lagrangian including fermion ghosts and auxiliary commuting fields, this could be an 
interesting topic for a separate publication. 



5 Conclusions 

In conclusion, we have seen that the leading order curvature terms as quadratic corrections to Einstein's 
gravity define a theory with a number of striking properties. The absence of van Dam-Veltman- 
Zakharov discontinuity in AdS space in the limit of the massive graviton is just a simple 

example shown above. Though one may have to check whether these results persist up to the loop level, 
such effects are much suppressed and would be insensitive to the present experimental observations, 
like gravitational bending. Nonetheless, one can expect vDVZ discontinuity in the graviton propagator 
to be absent even at loop level since the mass of the graviton is generated dynamically from the higher 
curvature terms. There is a clear message that a quadratic gravity itself deserves as a consistent theory 
of massive spin-2 field with sensible massless limit. This is the main thrust of this paper. We also 
obtained the general solutions to the linearized higher derivative field equations by specializing to the 
A = case in the Teyssandier gauge and explicitly identified the physical eight degrees of freedom 
of the theory. Finally we completed our discussion with a cursory inspection of the general gauge 
covariance by including the gauge fixing term. 

We argue that in any effective 4D theory obtained from the dimensional reduction of higher di- 
mensional theory both and Hq {Hq is the Hubble parameter) could depend on effective AD cos- 
mological constant. Since A ~ Hq, the ratio M^/A ~ M'^/Hq when remains finite, a massive graviton 
may contribute to gravity with a meaningful phenomenology. It is probable that a small fraction of 
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the gravitational interactions is associated with a massive graviton component, which could be even 
ghost-like, while still dominant is the component of massless gravitons, and this is indeed a general 
feature of four-derivatives gravitation. In brane world models, on the other hand, one may realize 
similar effects by having an ultra-light spin-2 particle with a very small coupling compared to massless 
graviton as pointed out in |8|, 29 1. 
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Appendix A: Linear Field Equations for i?^-Gravity 

For the mathematical clarity, we give the following variational formulae for different curvature and 
derivative terms. 

(i) Variation of G^i, 

6Gfj_i, = 5R^v — 2 dniySR — -jh^yR =^ SR^i, = 5G^v + ^ duu^R + ^ Vi^-^ 

(A.l) 

(ii) Variation of /^j. (= R [R^y - \ g^uR)) 

^I/iu = (^Rfiu — 2 9iJ.uR^ 5R + RSRpy — - hpyR? 

= Kg^, (y^h - V^V^hxp) + 2A (a^ + 2V(^V^/i,)A - VpV.hJ - AA^ V " \ 9puh^ 

(A.2) 

(iii) Variation of J^^u (= Rp.puaR'"^ - \ gp.vRpaR'"') 



6Jpu — SRppu "^R^ a + Rppva^R'"^ — t; 9puR'^'^6R, 



pu 



Rfipu Ra\ 2 9pvRp\Ru ^ h ^ ^ hpi,RpfjR^ 
I (a[') V + 2V(^V^/i,)A - V^V,/i) - 2A2 (^hf,, -^gp,h^ 
- 1 SRp, + I (v'/i - V^/iAp) ( A.3) 
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(iv) Variation of K^, (= R^p^sRu - \ Qi^vRpaixR""^^ - 2Rp ^R,x + 2R^x.pR^'') 
SKp, = 2i?(^ '"^i;,! 6R,)p^ ^ - i g^.uRP''^ JRp^x " + 2RP JR^p, ^ + 2Rpp,JRP'' - 4i?(^ P5R,)p 

(A.4) 



(v) Variation of VpV^R 

6 (VpV^R) = VpVJR - VxR 6Tp, ^ (A.5) 

where 

(vi) Variation of gp,i, V^R 

5 {gpuV^R) = V^Rhp, - gp^h^PVxVpR + g^, (vHR - gP'^STp^ ^ V xR) (A.7) 



(vii) Variation of V^i?^ 



pv 



S {V^Rpu) = V^6Rp, - hP'^VpVaRp, - gP'^dTp^ ^ VxRpu - {2VPR^,a6r,p " + Rp^dT.p + /i ^ z/) 

(A.8) 

Eq. (|2|) can also be written in the fohowing form (k^ = 1) 

(^Rpu - \gpuR^ + Agpu + 2aR {llp^ - -^g^uR^ + 2^ (^RppuaR"" - -^gpuRpaR^"^ 
+ 27 ^RR^u — 2Rppy„RP" + RppaxRu '"^^ — 2Rp ^Ryp — - gpuT^cB^ 

- (/3 + 2a) (V^V^ - gp^V^) R + pV^ (^Rp^ - ^ g^^R^ = , (A.9) 

where, as defined in the text, d = a — 7, /? = /? + 47, and the Gauss-Bonnet curvature squared 
^GB = R^ — + R^aXr ■ ^'^^ background solutions dH) , terms in the bracket multiplying with 
27 have a trivial contribution for the linear perturbations. While, on the Einstein backgrounds dH), 
the non-trivial contribution from the last two terms that involve covariant derivatives is 



'^(3 + 2aj {VpV^ - gp^V^) 6R + P V^6Gpu + 2A/3 V • (A.IO) 
Finally, using the results (i)-(iii), and ( A.10| ), we arrive to Eq. (11) given in the text. 
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